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ABSTRACT 

We compare the outer radius of the accretion disc in the intermediate-mass black hole 
candidate HLX-1 as estimated from the UV/optical continuum, with the values es- 
timated from its outburst decline timescales. We fit the Swift 2010 outburst decline 
lightcurve with an exponential decay, a knee and a linear decay. We find that the disk 
has an outer radius 10^^ cm < i?out ^ 10^^ cm, only an order of magnitude larger 
than typical accretion discs in the high/soft state of Galactic black holes. By contrast, 
the semimajor axis is « a few 10^"* cm. This discrepancy can be explained with a 
highly eccentric orbit. We estimate the tidal truncation radius and circularization ra- 
dius around the black hole at periastron, and impose that they are similar or smaller 
than the outer disk radius. We obtain that e > 0.95, that the radius of the donor star 
is < a few solar radii, and that the donor star is not at risk of tidal disruption. If the 
companion star fills its Roche lobe and impulsively transfers mass only around perias- 
tron, secular evolution of the orbit is expected to increase eccentricity and semimajor 
axis even further. We speculate that such extremely eccentric systems may have the 
same origin as the S stars in the Galactic centre. 

Key words: accretion, accretion discs - X-rays: individual: HLX-1 - black hole 
physics. 



. 1 INTRODUCTION 

'. The point-like X-ray source 2XMM J011028. 1-460421 

■ (henceforth, HLX-1 for simplicity) is the strongest 
, inter medi ate- mass black h ole (IMBH) candidate known to 

■ date IpWr cU ct al. 2009; Wiersema et al.ll2010l : iDavis et all 
• I2OIII : IServillat et all 120111 1. It is seen in the sky at a dis- 
, tance of « 8" from the nucleus of the SO galaxy ESO 243-49 

■ (redshift 2 — 0.0224, luminosity distance ~ 95 Mpc, dis- 
' tancc modulus ~ 34.89 mag; at this distance, 1" ~ 460 pc). 

, Its X -ray luminosit y and spectral varia bility (|Farrell et al.l 

■ I2OO9I : ICodet et all l2009i : IServillat et all I2OIII ) and its ra- 
dio flares detected in association with the X-ray outbursts 
l|Webb et al.ll2012l ) are consistent with the canonical state 
transitions and jet properties of an accreting BH. With a 
peak X-ray luminosity ~ 10*^ erg s~^, the BH mass required 
to be consistent with the Eddington limit is ^ 10^ M©. 
A similar value is obtained from spectral modelling of the 
thermal X-ray component, which is consisten t with emis- 
sion f rom an accretion dis c (|FarrelI et al ] |2009l : lDavis et all 
l201ll :[ Servillat et al.|[201ll l. If these BH mass estimates are 
correct, HLX-1 is way too massive to have been formed 
from any stellar evolution process. A more hkely scenario 
is that it is the nuclear BH (perhaps still surrounded by 
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its own nuclear star cluster) of a di srupted dwarf satel- 
lite galaxy, accreted by ESO 243-49 (|Mapelli et all I2OI2I : 
iKing fc DehnenI bOOSl ). HLX-1 has a point-like, blue opti- 
cal counterpart (_B ^ T/ ^ 24 mag near t he outburst peak; 
iFarrell et al]|2012l : [Soria et al.ll20ia . l2010l '). The presence of 
Hq emission at a redsh ift consistent with that of ESO 243-49 
(|Wiersema et al.ll2010l ) is perhaps the strongest argument 
for a true physical association. It is still debated whether 
the optical continuum emission is dominated by the outer 
regions of the B H accretion disc, or by a young star cluster 
around the BH (|Soria et al.ll2012l : iFarrell et alji2012l ). 

In the absence of phase-resolved dynamical measure- 
ments of the BH motion, we can use the Swift X-ray 
lightcurve properties to constrain the system parameters. 
The X-ray flux shows recurrent outbursts every m (366 ± 4) 
d (seen every late August in 2009, 2010, 2011 and 2012), 
due either to some kind of disc instability, or to a periodic 
enhancement of the accretion rate. Se veral alternat i ve sce - 
narios were considered and discussed bv lLasota et al.l l|201lf ) , 
who favoured a model in which enhanced mass transfer into 
a quasi-permanent accretion disc is triggered by the passage 
at periastron of an asymptotic giant branch (AGB) star on 
an eccentric orbit (e ^ 0.7). Since the publication of that 
work, the detection of the th ird and fourth consecutive out- 
bursts (see lGodet et al.||2012l for the first report of this year's 
outburst) has clinched the interpretation of the recurrence 
timescale as the binary period. Furthermore, additional op- 
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tical photometric results have been pubh shed, based on data 
from the Hubble Space Telescope (HST) (|Farrell et al.ll2012l ') 
and from the European Southern Observa tory (ESO)'s Very 
Large Telescope (VLT) llSoria et al.ll2012h. T hus, in this pa- 
per we revisit and update lLasota et al.l ( 201ll ')'s orbital mod- 
els and constraints in the light of the new results. 



2 SIZE OF THE ACCRETION DISC 

2.1 Predictions for a standard disc model 

At a distance of 95 Mpc, the characteristic size of the region 
responsible for most of the soft, thermal [kT ~ 0.2 keV) X- 
ray emission is ~ a few 10® cm (inferred from fits to XMM- 
Newton, Chandra, Swift spectra), and is consistent with 
being constant during the decline of in dividual outbursts , 
and over the three r e corded outbursts (|Farrell et al.l 
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Scrvillat ct al.' '20 1ll: iDavis et alj l201ll : ISoria et al.] 
2012 : FarrcU et al.[|2012h . This suggests that the soft X-ray 
emission traces the true inner radius of the disc, bounded 
by the innermost stable circular orbit around the BH. In- 
stead, much less is known about the outer disc radius, from 
UV/optical/IR observations; it is still debated how much 
of the blue optical emission comes from an irradiated disc, 
and how much from a possible cluster of young stars around 
the BH. If the disc is th e dominant UV/opt ical emitter, the 
HST a .nd VLT studies of lFarrell et"all l|2012l ') and lSoria et~al] 
(|2012l ) respectively, agree on an outer disc radius ~ 10 cm 
~ 1 arp. If a substantial contribution comes from unresolved 
young stars, we can take that value as an upper limit to the 
true disc size. A ratio of outer/inner disc radii ~ 10^ is sig- 
nificantly smaller than observed in transient Galactic BHs 
with Roche-lobe-filling donors, where typical outer radii are 
~ lO'^^ c m ~ a few 10^ times the innermost stable circu - 
lar orbit (|Zurita Heras et al.|[2oTll : iHvnes et~al]|2002l . lT99l ). 
This serves as a warning that we have to disentangle what 
scales with BH mass and what does not, when using scaled- 
up Galactic BH models to interpret HLX-1. While the inner 
disc depends directly on the BH mass, the outer disc de- 
pends mostly on the donor star and binary separation. 

There is an alternative way to estimate the outer disc 
siz e, based on the X-ray out burst decline timesc ale. FoUow- 
ing lKing fc Ritteij l{W9^ and lFrank et al.l (|2002l '). we assume 
that the outbursting disc is approximately in a steady state 
with surface density 



pH'. 



Mb 



(1) 



where Mbh is the central accretion rate and f the kinematic 
viscosity. When the whole disc from iiin to Tiout is in a hot, 
high-viscosity state, the total mass in the disc 

^ Note for the arXiv version: after this paper was accepted in 
MNRAS, a re-analysis of the HST data with a different disc irra- 
diation model, by Mapelli et al. (submitted to MNRAS), showed 
that they are consistent with an outer disk radius ~ 3.5 x 10^^ 
cm, for a viewing angle i = 45° ; the radius is smaller for a face-on 
view, which we consider a more likely scenario, given the narrow 
full-width-half-maximum of the Hq emission line (work by Hau 
&; Soria, in preparation). Even if we adopt this slightly larger up- 
per limit for the disc radius, the argument for a high eccentricity 
discussed in Section 3 remains unchanged. 
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where we have neglected other sources of mass loss from the 
disc apart from BH accretion. In Eq.(2), v is interpreted as 
an average value of the kinematic viscosity over the whole 
disc; in p ractice, we take the v alue of f near the outer edge of 
the disc (|King fc RitteJI 19981 ). Integrating Eq.(2), we obtain 
the well-known exponential decline for the disc mass 
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Mdisc = Mdisco exp(-3i^t/i?o 
and consequently also for the accretion rate 

Mbh = — ^2 exp(-3!/t/i?out). 

-"■out 



(3) 



(4) 



and the outburst luminosity L ~ Lx ~ 0.1i\fBHC . In sum- 
mary, we expect to see a luminosity 

Lx ~ Lx.o exp(— 3i^t/i?out), (5) 

where Lx,o is the value at the outburst peak, declining on a 
timescale 



-Rout /(Si'), 



(6) 



as long as the rate at which the disc mass is depleted during 
the outburst decline is much larger than any ongoing trans- 
fer of mass from the donor star. For the viscosity at the 
outer edge of the disc, we take the usual parameterization 
H (|Shakura fc Sunvaevll 19731 ) . where a is the viscos- 
ity coefficient in the hot state, Cs is the sound speed and H 
the vertical scaleheight. In the simplest, order-of-magnitude 
approximation, we can take an outer disc temperature « 10"* 
K (enough to keep it in the hot state), corresponding to 
Cs « 2 X 10® cm s~^, and a vertical height H « O.li?. This 
gives, from Eq.(6): 



Rou 



6 X 10 arc cm. 



(7) 



which we shall directly compare with the observations. 

If we adopt the Shakura-Sunyaev disc solution 
(|Shakura fc Sunvaevll973 : lFrank et al.ll2002l ). with Kramers 
opacity, we can write 



1.8 X 10"Q*/^Mf/^V-i/*<* cm^s-i 



1.4 X 10'" a 



16 4/5 



,>3/10 -l/4„3/4 2 -1 

m.3 cm S , 



where Mie is the accretion rate in units of 10^^ g s 



(8) 



mi 



is the BH mass in solar units, _Rio 
Then, from Eq.(6): 



Te ^ 5.2 X lO'^a'*^'' 
that is 

" (5:2^) 



M, 



4/5 



-3/10 1/4 p5/4 
2 //i.3 nj2 



16/25 ,V6/25^ 



i?out/(10'"cm), etc. 



(9) 
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(10) 



The peak luminosity Lx.o in Eq.(2) can be left as a 
purely observational parameter, or it can itself be expressed 
as a function of outer disc radius, viscosity, density and 
BH mass, if we assume t hat the outburs t is triggered via 
the dwarf- nova instability (|Gannizzdl 19931 ). (More precisely, 
if we assume that the outburst starts when the enhanced 
mass transfer due to periastron passage pushes the disk from 
the cold to the hot state.) In that case, the surface density 
at any radius immediately before the start of the outburst 
approaches the maximum value allowed by the S-curve in 
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the surface-temp erature phase space (|King fc RitteJ [19981 : 
iFrank et al]|2002l '): 



11 ,1 nl-05 -0.35 -0.86 -2 

11.4i?io m-^ g cm 

1 Q in2 nl. 05^-0. 35„-0. 86 „ „„-2 

i.^ X iu Ki2 TTig Q(. g cm , 



(11) 



where Oc ^ 0.01 is the viscosity parameter in the cold disc 
state. Taking for simphcity H ~ bR, where 6 is a constant 
~ 0.1, we can then express the maximum volume density at 
the start of the outburst as Pmax ~ ^ma^/H « T,^s.x/{bR), 
which is essentially independent of R, given the expression 
for Emax in Eq.(ll). It is then easy to integrate the total 
disc mass at the start of the outburst: 



2tt 



(2-Kb) 



RdR 

flout n3 

pR^dR^{2'Kb)^^ 



the disc mass at later times 



Mdis 



{2-Kb)pR:t 



exp(-3i^t/i?out), 



the accretion rate (cf. Eq.(4)) 

Mbh ~ {2%b){Ro^ty p) exp(-3!/f/i?out), 

and the peak luminosity 

Lx,o ~ {0.lc)(2TTb){R^^t^p) 

« (0.1c^)(27r6)i?out(aCsi/)(En.ax/-ff) 
^ 0.lQ:C^(27r6)i?outCsEmax, 



(12) 



(13) 



(14) 



(15) 



where Emax comes from Eq.(ll). 

So far, we have assumed that the whole disc is in the 
hot state; this is usually the case in the early part of an 
outburst, especially when the outer edge of the disc is kept 
in the hot state by X-ray irradiation. The exponential decay 
continues until the outer disc annuli can no longer be kept 
in the hot state, so that hydrogen recombines and viscos- 
ity drops. From that moment, the outer edge of the hot disc 
R\i < Rout - The contribution to the accretion rate and to the 
continuum X-ray/UV/optical emission from the outer (cold, 
low-viscosity) an nuli at Rh < R < Rou t becomes negligible. 
It was shown by iKing fc Ritteil l| 19981 ') that the central ac- 
cretion rate in this second phase of the decline is 



Mbh = MsHiti) [1 - C(t - ti)] , 



(16) 



where ti is the time after which the outer disc is no longer 
in the hot state, and C parameterizes the fraction of X-ray 
luminosity intercepted and thermalized in the outer disc. 
As C can be taken as a constant, Eq.(16) shows that the 
accretion rate and luminosity decline in the late part of the 
outburst is linear. Most importantly for our current purpose, 
the slope of the linear decline is such that 

^end — tl = Te (17) 

ijKing &i Ritteij|l998 l'). where tend is the (extrapolated) time 
in which the accretion rate and luminosity go to zero. 

Finally, we need to consider the case when there is on- 
going mass transfer M2 from the donor star during the out- 
burst. In that case, the asymptotic value of the luminosity 
in the exponential decline is not zero but L2 ~ 0.1(— M2)c^ 
(assuming a standard radiative efficiency « 0.1). Recalling 



that tl is the time when the lightcurve switches from an ex- 
ponential to a line ar decline, and defining Li = L{ti), Eq.(5) 
is modified as (Powell et al. I [2007. ) 

Lx = (ii - L2) exp(-3i/ {t - tl) + L2. (18) 



After the transition to a linear regime, the luminosity is 



1 



R^ 



■(i-ii) 



(19) 



Note that if —M2 > 0, the first derivative of t he luminosity 
is discontinuous aX t = ti |PoweireiaD[2003), because the 
gradient of the exponential decay is 



Lx{ti) 



3u 

r>2 



Li 1 + 



0. 1A/2C^ 



that is flatter than the gradient of the linear decay 
3u 



Lxiti) = 



■Li. 



(20) 



(21) 



Therefore, the exponential-to-linear transition is often re- 
ferred to as the "knee" in the lightcurve of transient X-ray 
binaries. 



2.2 Comparison with the observations 

We shall now fit the X-ray lightcurve to obtain two inde- 
pendent estimates of the viscous timescale To, from the ex- 
ponential and the linear regime, and use them to constrain 
Rout from Eq.(lO), or, using a simpler approximation for 
the scale-height, from Eq.(7). We shall then derive an inde- 
pendent estimate of Rout from the expression for the peak 
luminosity in Eq.(15). We st udied the publicly a vailabljfl 
Swtft X-Ray Telescope (XRT; iBurrows et"aLll2005h data for 
the 2010 outburst, because it is the same outburst for which 
we o btained constraints on -Rout from t he optical contin- 
uum l|Soria et al.ll2012l : iFarrell et al.ll2012l'). We used the on- 
line S wift/XRT data product generatoiFI IjEvans et al.ll2007l . 
I2OO9I ) to extract a lightcurve in the 0.3-10 keV band. We fit- 
ted the lightcurve with an initial exponential decay, a knee 
and a linear decay (Fig. 1). The shape of the X-ray out- 
burst lightcurve of HLX-1 is remarkably similar to those of 
several transient G alactic X-ray binari es (BHs and neutron 
stars), modelled by I Powell et al.l l|2007h . which were success- 
fully used to constrain the size of their accretion discs. 

For the exponential part, we obtain a best-fitting 
timescale To = RLt/{3iy) = 3.7^1° x 10*^ s (90% confidence 
limit). For the linear part, we have Tc = 3.5lQg x 10® s. 
Assuming a peak luminosity ~ lO'*^ erg s~^ in the 0.3-10 
keV band, and a bolometric luminosity a factor of 2 higher, 
implies an accretion rate M « 2 x 10^^ g at stan- 
dard efficiency. The viscosity par ameter a < 1, and more 
likely a ~ 0.3 (|Frank et al.l 120021 ). From Eq.(lO), this im- 
plies i?out ~ lO'^^ cm, only very weakly dependent on BH 
mass and accretion rate. Using the approximation in Eq.(7), 
we also obtain -Rout ~ 10^^ cm. From the peak luminosity 
(Eq.(15)), for a ~ 0.3 we obtain Rout ~ 4 x 10^^ cm. 

We also analyzed the lightcurves for the 2009 and 2011 



^ http:/ /www. swift. ac.uk/user_objects 

Including the new treatment of the vignetting correction, in- 
troduced after 2011 August 5. 
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Figure 1. Swift/XKT lightcurve of the 2010 outburst, fitted with 
a standard X-ray transient model (exponential decay, knee, linear 
decay) . 



outbursts (the 2012 outburst was still ongoing as this pa- 
per went to press). They are more noisy, less easy to inter- 
pret in terms of exponential and linear branches. However, 
for both of them it is possible to estimate an e-folding de- 
cline timescale, roughly corresponding to the exponential 
timescale determined for the 2010 outburst. The timescales 
are 5 X 10® s and s:; 3 x 10** s for 2009 and 2011, re- 
spectively, and the peak luminosities are approximately the 
same in all three outbursts. Thus, we also estimate an outer 
radius ~ 10^^ cm in the 2009 and 2011 outbursts, in the 
standard disc approximation. 

Those values are almost one order of magnitude smaller 
than what we estimated by assuming that most of the con- 
tinuum emission comes from the hot disc; and the latter was 
already a surprisingly small radius compared with the binary 
system parameters (see Sect ion 3 ). Wc^notc that bot h the 
HST and VLT observations ijSoria et aT.ll2012l : iFarrell et al.l 
l2G12h were taken during the exponential part of the decline, 
that is when the whole disc was in a hot state. Therefore, the 
estimates of -Rout from the optical/UV continuum should be 
comparable to those from the X-ray lightcurve. New optical 
observations early in the next outburst will hopefully allow 
us to measure the true outer disc size and luminosity. 



3 ORBITAL PARAMETERS 

We can now compare the size of the disc estimated from 
optical and X-ray flux measurements (10^^ cm < iiout < 
10^^ cm) with the characteristic size of the binary system. 
Because of the sharpness of the outbursts rise and decline, 
reminiscent of Galactic X-ray binaries, we assume that the 
BH is accreting from a single donor star rather than AGN- 
like gas inflows. We also assume that the outburst recurrence 
timescale ~ 370 d corresponds to the binary period. Then, 
the semimajor axis a of the binary is 



a = 1.50 X lO^^m^''^ 



+ 



1/3 p2/3 



Pyr ' cm 



(22) 



(jNewtonllTeSTl ). where q = M2/MBH and m = (A/bh + 
M2)/Mq = M/Mq. Typical values for HLX-1 in the 
intermediate-mass BH scenario are q ~ 10~^ and m^^^ ~ 



10-20. Therefore, in the most accepted scenario, the semi- 
major axis is at least 10, and possibly up to 100 times larger 
than the disc radius. This misma tch clearly suggests an ec- 
centric orbit l|Lasota et al.ll201ll ). in which the character- 
istic disc size is determined by the periastron separation 
Rpei = (1 — e)a, with eccentricity e > 0.9. 

The amount of mass transferred to the BH in each 
outburst suggests that the donor star overflows its instan- 
taneous Roche lobe every time it passes at periastron. In 
Roche-lobe mass transfer systems, the outer edge of the 
accretion disc is generally identified with the largest sta- 
ble non-intersecting orbit (tidal truncatio n radius Rt)- For 
mass ratios M2/M1 < 1 Rt 0.48a JPaczvnsk: Il977l : 
iPapaloizou fc PringldlmTl : I WhitehurstI 19881 : 1 Warnedl 19951 ) . 
If unstable orbits are also allowed, the disc may expand up to 
-Rt ~ 0.60a/(l + g) (Warner 1995, and references therein). If 
we take the periastron distance as the instantaneous binary 
separation during the phase of Roche-lobe overfiow, this cor- 
responds to an expected disc size -Rout ~ 0.6(1 — e)a. For an 
observed disc size -Rout ~ 10^'^ cm, the tidal radius condi- 
tion would require an eccentricity e « 0.89 for a BH mass 
= IOOOM0, and e ^ 0.95 for a BH mass = IO^Mq. If we 
take the lower bound to our observed disc size -Rout = 10^^ 
cm, we need e « 0.989 or e ~ 0.995, respectively. 

We can argue that the tidal truncation constraint is not 
relevant to the case of HLX-1, where mass transfer may oc- 
cur impulsively near periastron, and the timescale for the 
disc to expand to its tidal truncation radius is similar to 
the timescale for the disc matter to be accreted and for the 
binary orbit to expand after periastron. In other words, in 
HLX-1 the disc may look small because it did not have time 
to grow to its tidal truncation radius. Instead, the circu- 
larization radius provides a stronger lower limit to the pre- 
dicted disc size, and is applicable to any system where mass 
transfer occurs through the Lagrangian point Li. 

The circularization radius -Rcir is defined via the con- 
servation of angular momentum equation 



Vtf, {Rcir) Rcir — (^LlPpcr)^ ^ (Rpcr) ■ 



(23) 



where is the orbital velocity of the accretion stream 
around the BH, J^Li-Rpcr the distance between the BH and 
the Lagrangian point Li, and Q{Rpcr) the angular veloc- 
ity of the donor star at periastron. Here, we must be care- 
ful not to use th e well-known fitti ng formula for Xli ~ 
0.500-0.227 log g (jPrank et al.ll2002l ) because it applies only 
for g > 0.1 and for circular orbits. Instead, wc need to com- 
pute Xli from Eq.(A13) of ISepinskv et alj (,2007 ). valid for 
eccentric orbits: 

1 .2 -^-/(l-^Li)(l + g)(l + e) + 1^0.(24) 

(1— Ali) A-Li 

Here, / is the ratio of the rotational angular velocity of the 
donor star to its orbital angular velocity, at periastron. It 
parameterizes the degree of tidal locking; a star with / = 1 
is rotating synchronously with the orbit, at periastron. A 
list of Ali solutions for characteristic values of q and e is 
given in Table 1. The orbital velocity of the accretion stream 
around the BH 



( GMbh \ 

\ Rcir J 



1/2 



(25) 



and the orbital angular velocity at periastron 
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f=0 









e = 0.80 


e = 0.90 e = 0.95 


e = 0.99 


1 = 


10- 


1 






0.695 




Q = 


10- 


2 






0.843 




q = 


10- 


3 






0.924 




1 = 


10- 


4 






0.964 




1=1 


1 = 


10- 


1 


0.732 


0.734 


0.735 


0.735 


q = 


10- 


2 


0.868 


0.869 


0.869 


0.869 


1 = 


10- 


3 


0.937 


0.938 


0.938 


0.938 


1 = 


10- 


3 


0.970 


0.971 


0.971 


0.971 



Table 1. Distance Xli between the BH and the Li point (as a 
fraction of the periastron distance ijpcr), in the parameter range 
of interest for HLX-1; from Eq.(24). 



(l + e 



,1/2 



(1 _e)3/2 



G(A/bh + M2) 



1/2 



(26) 



For the sake of our numerical estimate, we take Xli = 0.95, 
a typical value in the range of parameters thought to be rel- 
evant for HLX-1 (Table 1). Then, substituting into Eq.(23), 
we have: 



(0.95fipe0^ jl + e) G (Mbh + M2 

GA/bh (1 - e)3 
0.81(1 -e^)(l + g)a 
1.2 X lO^mi/^'d 



(27) 

Assuming that the BH mass is in the range ~ lO'^-lO'^M©, 
Eq.(27) gives us a strong constraint on e, by imposing that 
the circularization radius is smaller than the observed outer 
disc size. For example, assuming Mbh = 5 x 10'^ Mq, a cir- 
cularization radius _Rcir = 10^'' cm (at the upper end of our 
disc size estimates) requires e ~ 0.97; for Rdi = W^^ cm (at 
the lower end of our disc size estimates), e « 0.997. The cor- 
responding periastron distance^ -Rpcr = (1 — e)a ^ 6 x 10^^ 
cm for e = 0.97, and Rper « 6 x lO" cm for e = 0.997. 

The distance between Li and the centre of the donor 
star places an upper limit on its radius; characteristic values 
of its instantaneous volume- av eraged Roche lo be radius at 
periastron can be obtained from lEggletonI (|l983h . For typical 
values q ~ 10~*-10~^, the secondary gets squeezed to a 
radius « (0.02-0.05)i?per. Since we have assumed that the 
secondary fills the Roche lobe and dumps mass into the BH 
only near periastron, this radius must also be similar to the 
size of the donor star. For example, if g = 2 x 10~* and e = 
0.97, the donor star must have a radius ~ 'IRq, consistent 
with main sequence and subgiant stars. 



* Note that for e > 0.2, the periastrion distance between the two 
stars is always smaller than the circularization radius around the 
accreting primary. This is because the angular momentum of the 
secondary at periastron, and of the matter transiting through the 
Li point, is larger than the angular momentum of a circular orbit 
with the same binary separation. It is not a problem, because by 
the time the accretion stream has completed a full orbit around 
the BH and formed a ring, the secondary has moved away from 
periastron and the primary's Roche lobe has widened. 



The values of e estimated here from disc size argu- 
ments are much mo re extreme than what was suggested in 
iLasota et all ((201J). They may seem implausible, knowing 
that tidal forces tend to c i rcular i ze orb its in X-ray binaries. 
However, ISepinskv et all (|2007l . \200^ ) showed that in the 
case of a donor star that transfers mass impulsively only at 
periastron, with g < 1 — 0.4e + 0.18e^, the secular evolution 
of the orbit leads to an increase of both eccentricity and 
semimajor axis, even when the opposite effect of tidal forces 
is taken into account. If tidal circularization is neglected, the 
eccentricity increases as 



{e} = -^{l-er^'il-e)iq-l) (28) 
and the semimajor axis (and hence the binary period) as 

2')'/'(9-l). (29) 



TT M2 



It is then easy to show from Eq. (28,29) that the periastron 
distance Rpcr = (1 — e)a, and therefore also the size of the 
secondary's Roche lobe at periastron, remains unchanged. 

Finally, we need to assess what donor stars can sur- 
vive on such eccentric orbits with small periastron distance, 
avoiding tidal disruption. The condition for survival is that 
the periastron distance Rpcr is larger than the tidal disrup- 
tion radius ()ReesT988. ) 



i?td ~ 5 x lO^mg/^ {R2/Rq) {M2/Mq)-^^'^ 



(30) 



By substituting R2 < 0.05i?por into Eq.(31), we can re-cast 
the tidal survival condition as 

M2 > 4.6 X 10~^Mbh, (31) 
easily satisfied in the likely mass range of HLX-1. 

4 CONCLUSIONS 

We compared the estimates of the disc size from the opti- 
cal continuum flux {Rout ^ 10^^ cm) with those obtained 
by fitting the X-ray luminosity decline after an outburst. 
We found that, at least for the 2010 outburst, the de- 
cline displays the standard sequence of exponential phase, 
knee, linear phase often seen in Galactic X-ray binaries; this 
strengthens the interpretation that the thermal X-ray emis- 
sion in HLX-1 comes from a disc, and the decline timescale 
corresponds to its viscous timescale. For all three outbursts 
observed to-date, the timescale is very short, ~ 4-8 weeks. 
This is similar or only slightly longer than what is typi- 
cally observed in Galactic X-ray transients, despite that fact 
that both the orbital period and the BH mass (and, hence, 
the semimajor axis) of HLX-1 are > 100 times larger. The 
outer disc radius estimated from the viscous timescale is 
Rout ~ lO'^'^ cm, if the viscosity parameter is similar to the 
values usually estimated for Galactic BH transients in a high 
state. We cannot rule out that the fast accretion of the disc 
matter in HLX-1 may be partly due to an effective viscosity 
Q^cff 1, higher than in the Shakura-Sunyaev prescription. 
But we argue that even if we assume the upper disc size es- 
timate -Rout ^ lO'^^ cm, a highly eccentric orbit is required 
to explain the small disc size. 

To quantify the eccentricity, we calculated the charac- 
teristic length-scales of the binary system, as a function of 
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BH mass and eccentricity. If the disc extends at least as 
far as the circularization radius (as is usually the case in 
X-ray binaries with Roche- lobe mass transfer), we obtain 
that _Rcir ~ (1 — e^)a, and therefore e > 0.95 for a BH 
mass > 10'^ Mq. We argued that X-ray binaries with such 
extreme values of e are the most likely evolutionary endpoint 
of systems with g <C 1 and a moderately eccentric initial or- 
bit, such that Roche- lobe-overflow mass transfer occurs only 
impulsively near periastron. Secular evolution will tend to 
make the orbit more and more eccentric, by increasing the 
semimajor axis and the binary period, at constant perias- 
tron distance. 

The small periastron distance required to explain the 
HLX-1 observations sets an upper limit to the current ra- 
dius of the donor star -R < a few Rq, ruling out supergiants, 
red giants and AGB stars. Possible donors are main sequence 
(B type or later) or subgiants. The compactness of the donor 
star, and the fact that secular orbital evolution due to mass 
transfer will not change the periastron distance, imply that 
the companion star in HLX-1 is not at immediate risk of 
tidal disruption, and will not be in the near future. In other 
words, we are not observing HLX-1 in a peculiar moment 
of its evolution, immediately prior to tidal break-up of the 
donor star. HLX-1 appears to be a stable system, with a life- 
time for X-ray outbursts determined primarily by the mass 
transfer timescale from the donor; at a rate ~ 10~^A{q yT~^ 
(averaged over the binary period), it may last for another 
~ 10^-10^ yr, during which its semimajor axis and binary 
period (and hence, interval between outbursts) will continue 
to increase. 

Eccentricities > 0.95 may seem implausibly extreme, 
but there is at least one class of stellar objects where they 
are the norm: S stars observed on highly eccentric orbits 
within 0.01 pc of the Galactic nuclear BH (Alexander 200^; 
iGillessen et"al]|2009f ). A possible scenario for the origin of 
Galactic S stars is the tidal disruption of a stellar binary 
system near the BH, which produces an escaping, hyper- 
velocity star, and a more tightly bound star on a v ery eccen- 
tric o rbit, theoretically as high as e ~ 0.99 (.Lockmann et all 
|2008| ). Observationally, the most eccentric, bound S star 
for which orbital par ameters have bee n reliably determined 
has e ~ 0.96 (Gille ^n et al.l |2009| ). We speculate that 
intermediate-mass BHs in star clusters may also capture 
stellar companions on very eccentric orbits through a similar 
process. 
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